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Abstract 

We study exclusive radiative decays of the T using soft-collinear effective theory and non- 
relativistic QCD. In contrast to inclusive radiative decays at the endpoint we find that color-octet 
contributions are power suppressed in exclusive decays, and can safely be neglected, greatly sim- 
plifying the analysis. We determine the complete set of Lorentz structures that can appear in the 
SCET Wilson coefficients and match onto them using results from a previous calculation. We run 
these coefficients from the scale Mx to the scale A ~ 1 GeV, thereby summing large logarithms. 
Finally we use our results to predict the ratio of branching fractions B{T ^ f2) / B{J /ip 7/2)) 
B^J/ip — > 7/2)/-B(V'' — *■ 7/2)) and the partial rate for T — > jnir. 
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I. INTRODUCTION 



In a recent series of papers the differential decay rate for the decay T — 7X has been 
studied in the "endpoint" region where the decay products have a large total energy of order 
the T mass (Mx), and a small total invariant mass squared of order AMx, where A ~ 1 GeV 
is the typical hadronic scale 0, S M S • An important tool in this analysis is the soft- 
collinear effective theory (SCET) p, H Hl^! which is a systematic treatment of the high 
energy limit of QCD in the framework of effective field theory. Specifically SCET is used 
to describe the highly energetic decay products in the endpoint region. The heavy b and b 
quarks which form the T are described by non- relativistic QCD (NRQCD) HI El. 



The soft-collinear effective theory is not limited to applications involving inclusive pro- 
cesses. In fact SCET has been extensively applied to exclusive decays of B mesons into light 
mesons BlIlliiHHIillillllSEllH Here we use similar techniques to study 
exclusive radiative decays of the T. We make use of some of the results derived in the anal- 
ysis of inclusive radiative decays in the endpoint region, but the analysis of exclusive decays 
is complicated by the existence of two different collinear scales. This necessitates a two-step 
matching procedure [T3|. In the first step one matches onto SCETj which describes collinear 
degrees of freedom with typical offshellness of order a/AMx, as is appropriate for inclusive 
decays in the endpoint region as discussed above. In the second step SCETi is matched onto 
SCETn, which is appropriate for exclusive processes since it describes collinear degrees of 
freedom with typical offshellness of order A. 

The analysis of T decay is further complicated by the existence of two types of currents: 
those where the bb is in a color-singlet configuration and those where it is in a color-octet 
configuration. The octet operators are higher-order in the combined NRQCD and SCET 
power counting, so one might suppose that t hey can b e dropped. However, the octet currents 
have a Wilson coefficient which is order a/ asiM-^) while the singlet current has a Wilson 
coefficient of order as{Mr)- The additional suppression of the singlet Wilson coefficient is 
enough so that both color-octet and color-singlet operators must be included as contributions 
to the inclusive radiative decay rate in the endpoint region j3, |^ . 

In this work we show that in exclusive decays the octet currents are truly suppressed 
relative to the singlet current and can be neglected. We then determine the minimal set of 
color-singlet currents which can arise and fix their matching coefficients in SCETj. We run 
this current to the intermediate collinear scale fic ~ \/AMr and match onto SCETn. Our 
expression for the decay rate agrees to leading order in the twist expansion used in previous 
work in QCD [i^ |2^. Finally we use our results to make a prediction for the ratio of 
branching fractions B{T — >• •y / B {J / ip — > 7/2), B{J/ip — > '~^f2)/B{ip' — > 7/2), and analyze 
the decay T 77r7r in the kinematic regime where the pions are collinear. 



II. POWER COUNTING 



A. Inclusive Decays 

The first step is to match the QCD amplitude for a hb pair in a given color and spin 
configuration to decay to a photon and light particles onto combined SCETi and NRQCD 
currents. The SCET power-counting is in the parameter A ~ a/A/M, where M = 2mfe, while 
the NRQCD power-counting is in f , the relative velocity of the bb pair in the T. Numerically, 
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FIG. 1: Matching onto operators in the effective field theory with one and two gluons in the final 
state. The currents on the left have a color-octet bb in either a ^Sq or ^Pj configuration. The 
matching for a color-singlet bb pair in a ^Si configuration is shown on the right. 

A ~ f ~ 1/3. The matching is shown graphically in Fig. ^ The effective theory operators 
can be classified into those with the bb in a color-octet configuration (shown on the left-hand 
side of Fig. P) and those with the bb in a color-singlet configuration (shown on the right-hand 
side of Fig. H}. The leading octet operators can be further subdivided into a those where 
the bb is in a ^5*0 configuration and those where the bb is in a ^Pj configuration. The octet 
^5*0 operators are 

U8,'So) = Y,Cf"'%M,^^)^^^xlpBl^p, (1) 

i 

where 

Bl = --W^[P'^+gM'n,,)±W. (2) 
9s 

The operator projects out the label momenta in the perpendicular direction The 
sum in Eq. (jT)) is over all possible Lorentz structures denoted by F^^, and Cf' ^°{M,iJ,) is 
the corresponding matching coefficient for each structure. The octet ^Pj operators are 

U8,'Pj) = ^Cp-'(M,/i)Fj,^^,xtp52A-rA-'T^V^p, (3) 

i 

where A is a Lorentz boost matrix. Each of these color-octet operators scales as 0{X) in 
SCET. The NRQCD power-counting has the ^5*0 octet operators scaling as 0{v^); however, 
this operator has an overlap with the T state beginning at (9(f^). Thus the ^5*0 operator 
contributes at order v^X to the T radiative decay rate. The ^Pj octet operator has NRQCD 
scaling 0{v'^), but overlaps with the T at order v. Thus the total power-counting of the 
^Pj contribution is 0{v^X), which is the same as the ^5*0 octet operators. The leading order 
matching coefficients for both are 0{^/aJ^W)). 
The color-singlet operators are 

J,a, 'Si) = J2 r^xtpA ■ a'i;^TT{Bl C^'''^\m, V+) 3^] , (4) 

i 

where V+ = "P^ + "P, with V = n ■ V . These operators scale as 0{X^) in SCET and 0{v^) 
in NRQCD. The leading matching coefficients are 0{as{M)). Thus the ratio of color-octet 
to color-singlet contributions in inclusive radiative T decay scales as: 

octet V . . 

~ ~ . (5) 

singlet X^/ajM) ^as{M) 
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FIG. 2: An example Feynman diagram of a time-ordered product in SCETi that matches onto an 
operator in SCETn that has a non-zero overlap with the final state collinear meson. 



B. Exclusive decays 

The situation changes when one considers exclusive decays. The currents we just discussed 
are SCETj currents where the typical invariant mass of the collinear degrees of freedom is of 
order /i^ = V MA. In order to have overlap with the meson state we must match onto SCETn 
currents where the typical invariant mass of collinear particles is 0{A). Furthermore, the 
interpolating field which annihilates the meson state in SCETn is defined to consist only of 
collinear fields in a color-singlet configuration jl^. Given these considerations it is simple 
to match the color-singlet operator in SCETj to an operator of identical form in SCETn. 
However, the matching of the octet contributions from SCETi onto SCETn is more involved. 

Before we consider the matching of the octet contributions from SCETj onto SCETn 
we turn our attention to the scaling of these contributions in SCETj. In order to produce 
a final state consisting only of collinear fields in a color-singlet configuration we need an 
interaction which changes the ultrasoft (usoft) gluon into a collinear gluon (as shown on the 
left-hand side of Fig. This term in the SCET Lagrangian is power suppressed by A so 
that the time-ordered product of the octet current with the collinear-coUinear-usoft vertex 
scales as (9(A^) in the SCETj power counting. In addition, the exchanged gluon introduces 
an extra factor of the coupling constant at the matching scale: as{fic)- Including these 
factors, the time-o rdered product of octet currents with the subleading Lagrangian scales as 
0{as{i^c)^/c(s{M)X'^v^), and the ratio of time-ordered products to the singlet contribution 
is 

^-^^^-0.05, (6) 
singlet ^a,(M) 

for the bottomonium system. For charmonium the above ratio is about 0.2. This result is 
very different from the result for the inclusive decay given in Eq. (jH)). In SCETj the octet 
contribution to exclusive radiative T decay is not only suppressed in the limit v,X — >■ 0, but 
numerically suppressed by a factor of ~ 10 for typical values of the parameters. This is the 
same order of suppression we expect from higher order SCET and NRQCD corrections; thus, 
we should be able to safely neglect the color-octet contribution in SCETj. However, before 
we can neglect the octet contribution in our analysis we must show that the suppression of 
the octet piece holds after matching onto SCETn. 

We first turn our attention to the simpler calculation: matching the color-singlet operator. 
In SCETi we perform the field redefinition 

A„-.rAWFt, (7) 
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which decouples usoft from colhnear degrees of freedom. Under this field redefinition 

Bi ^ y^f , (8) 

and the F's cancel in the trace of the color singlet operator given in Eq. (@]). Thus we 
match the SCETj operator after the field redefinition onto an operator in SCETn of a form 
identical to that in Eq. 



where /ic = V MA is the SCETi-SCETn matching scale, and the subscripts indicate SCETn 
fields. From now on we drop the subscripts. In SCETn the power-counting parameter is 
7] ~ A^, and the SCETn color-singlet operator in Eq. Q is O^v^r]"^). The short- distance 
coefficient is inherited from SCETj and is 0{as{M)). 

The matching of the color-octet current is more complicated. In order to match onto 
an SCETn operator with color-singlet collinear degrees of freedom we must consider time- 
ordered products where a usoft gluon radiated from the bb pair is turned into a final state 
collinear degree of freedom. An example of such a diagram is given in Fig. |21 Two collinear 
gluons are required for the collinear final state to be color-singlet. One of the collinear 
gluons comes from the octet current, and the other can be produced by pulling a gluon out 
of the bbg Fock state of the T, and kicking it with a collinear gluon from the current. This 
requires a colhnear-coUinear-ultrasoft coupling which first appears at order A in the SCETi 
Lagrangian [HHH: 



= ^ Tr { [^V^, zD^n [tV,, zDi] } , (10) 

where = + n- Dugfi^^ /2. The decay amplitude comes from a time-ordered product 
of the color-octet current and £^cg , or a time-ordered product of the color-octet current, 
C^cg , and a leading order gluon interaction. Though our result will hold for either type of 
time-ordered product we will, for the sake of concreteness, only consider the former: 

Ts = J rf^a;T{j(8,-)(0),/:W(x)}, (11) 

where the dot stands for ^5*0 or ^Pj. In the time-ordered product two gluon fields are 
contracted to form the internal propagator in Fig. |21 which scales as 1/A^. We require 
two uncontracted fields (in a color-singlet configuration) so that we can match onto an 
SCETn operator in the form Tr[y4;!-jy4^] which annihilates the final state collinear meson. 



In this example the leading contribution is an n ■ gluon field in one of the Wilson lines in 
J(8, ■) contracted with axin- An field in C^J . After the contraction what remains in £^cg is 

j,Tr{[gT^,A^:][P,Ai]}, (12) 

which scales as (A^)(1)(A) = A^. Note we now have the correct field content for the operator 
shown on the right-hand side of Fig. |21 there are two outgoing A']_ fields, one from C^cg and 
one from J(8, •), in a color-singlet configuration, and an incoming soft gluon field also from 

t~-cg ■ 
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Next we decouple collinear and usoft in SCETj through the field redefinition in Eq. ((Tj). 
This introduces factors of Y and into our expressions. When matching onto SCETn these 
become soft Wilson lines S and S'^. Since these Wilson lines do not affect the power counting 
we ignore them. Now we can match onto a convolution of SCETn operators with SCETj- 
SCETii matching coefficients. Since these arise from integrating out the internal collinear 
propagators they scale as for each propagator. In our example there is one propagator 
so the matching coefficient scales as (9(A~^), which is 0{ri~^) in the SCETn power counting 
(remember rj ~ A^). Since the SCETn operator has two Ajr^ fields each scaling as r], and a 
soft field scaling as 1], it scales as 1]^. Combining the scaling of the SCETn operator with 
the scaling of the SCETj-SCETn matching coefficient gives an 0{as{fJ'c)i]'^) contribution. If 
we include the order NRQCD scaling from the heavy sector, and the 0{^^/a^{M)) contri- 
bution from the QCD-SCETj matching coefficient the color-octet contribution to exclusive 
decays scales as 0{v^7]'^ as{M)as{fic)) in SCETn. Taking the ratio of the color-octet to 
the color-singlet contribution to exclusive T decay in SCETn we find: 

octet v^asific) 

~. — ; — ~ — , , (io) 

smglet ^as{M) 

which is the same scaling we found in SCETi. Thus we can safely neglect the color-octet 
contributions at this order. 



III. COMPLETE BASIS OF COLOR-SINGLET MATCHING COEFFICIENTS 

Now that the color-octet contribution has been eliminated we determine a complete basis 
of Lorentz structures T'^^p^iv ^^^^ can appear in the color-singlet matching coefficient in 
Eq. (0)). At leading order in as{M) only one Lorentz structure was found to be non-zero 0: 

CI (M, uj)T^ps^ = -^j^9a(39,^s . (14) 

However, at higher order other Lorentz structures may appear. These coefficients can be 
constructed from the set: 

{gi„y,nf,,nf,,Vf,} , (15) 

where v is the four- velocity of the T, under the restriction that F^^^^ satisfies the appropriate 
symmetries. For example, the full theory amplitude is parity even, as is the effective theory 
operator, meaning that the matching coefficient must also be parity even. As a result the 
epsilon tensor is not included in Eq. (jl5|) . 

We treat v as an object independent of n,n ji^, and use n"^ = = 0, n ■ n = 2, and 
= 1. Before we write down all the possible operators which can appear we note some 
simple properties that will make our task more manageable. First we note that Taps^i must 
be symmetric in a and (3. To see this consider the object: 

J2ci''''\M,u;)Tl,,^Tr{Bl5^^p^B'j . (16) 

UJ 

which is the collinear part of the color-singlet operator where a sum over u has been in- 
troduced. First interchange a and (3, and then use the cyclic nature of the trace to switch 
the two B± fields. Note, however, there is a projection on these fields from the operator 
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in the Kronecker delta involving Vj^. Since V+ = + V this operator projects out mi- 
nus the label on 5° and projects out the label on To preserve this relationship 
when the order of the fields is switched we must let S^^p^ ^uj,-p+- By letting u —uj 
we have 5_^_p^ = S^^p^, and the operator goes into itself. However, the Wilson coeffi- 
cient is now C^^' ^^\m, —u). To demonstrate that Eq. (fTBj) is symmetric under a ^ {3 

we must show that C\^' ^^\m,uj) is even in uj. We use charge conjugation for this. The 
heavy quark sector of the operator has charge conjugation C = — 1 as does the photon. As 
noted in Ref. 0| two gluons in a color-singlet configuration must have C even. Since QCD 
is charge conjugation conserving the product of operator and coefficient in Eq. (fTH|) must 

also be C even. This is the case if the matching coefficient C\^' '^^■'(M, a;)r^^^^ is C even. 
Following Ref. js^, under charge conjugation the above product of operator and coefficient 
goes to itself with uo — > —uo in the coefficient function. Thus charge conjugation implies 

as a result Eq. is symmetric in a and (3. Second, 
any vs appearing in F^^^^ gives zero contribution to the operator, since v ■ A = 0. Third, 
UayUa (and by symmetry njsjfip) appearing in the operator also gives a zero contribution 
since these indices contract with indices on the B± field. Finally, we use reparameterization 
invariance (RPI) of SCET jiil. l3ll|. The terms satisfying these requirements are 



+ C4 



9ci, \9fi& + 9 Pp. 



n ■ V 



So far we have allowed v to be an arbitrary vector. Now we restrict ourselves to a frame 
where = 0. Furthermore we are interested in the case where the photon is real, so 
we can restrict the photon to have transverse polarizations. This leaves only two linearly 
independent terms: 

^ Cf ''^'^ (M, uj)r^p^^ = a\9ip9i^ + af {9is9pp + 9ip9ps - 9^(59^) , (18) 

i 

where 9^ = 9^^ — {n^n'^ + n'^n'^)/2. The first term projects out the trace part of the 
Tr(i?"5^ p^i?^) operator, while the second term projects out the symmetric traceless com- 
ponent. Since the Lorentz symmetry in the perpendicular components is not broken in 
SCET these two terms do not mix under renormalization. The leading-order matching fixes 
the coefficients ai and 02 to order Ug'- 

a?(P+;/i = M) = ^|^, af(P+;/. = M)=0. (19) 

Since there is no mixing, 02 = + 0{a1{M)) at all scales. Note this matching assumes the 
T states are non-relativistically normalized: (T(P')|T(P)) = 5^{P — P'). 

In addition to gluon operators we must consider the basis of all possible quark operators 
which can appear in radiative T decays 

J^(l, ^Si) = J2 xtpA ■ a'^^^„,p,Wr^siM, V+)W^^n,p, . (20) 
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The basis of Dirac structures, '^75, -^7^}, was given in Ref. S^l, and the most general 
basis of quark operators can then be constructed out of these Dirac structures and the set 
{^afBfiu, gfii',nfj,,n^,Vf^}. Using the symmetries of SCET and RPI we find 

E = 4^9U + «^|75i + ai|7>. • (21) 

i 

The first term transforms as a scalar, the second term transforms as a pseudoscalar, and the 
third as a vector. The matching coefficients at the scale [i = M for the quark operators in 
radiative T decay are all zero at leading order in perturbation theory, but the scalar quark 
operator mixes with the scalar gluon operator through renormalization group running and 
can be generated in this manner. The pseudoscalar and vector term do not mix with the 
scalar gluon operator due to Lorentz symmetry and will not be generated at this order in 
the perturbative matching. 



IV. DECAY RATES & PHENOMENOLOGY 

We now consider the phenomenological implications of our analysis for exclusive radiative 
decays of quarkonium into either a single meson or a pair of mesons which are collinear. The 
(n + l)-body decay rate is given by: 

r(T-.7F„) = ^ f -^-^f]^-^i27T)^6\P-q-yp,) 

x|(7(g)pi...p„|JM,|T(P))p (22) 

where J is the QCD current, A is the photon field, and F„ denotes an exclusive final 
state consisting of n collinear particles. We consider only decay rates where the final state 
momenta Pi are all collinear with combined invariant mass m\ = {YllVi)"^ ~ ^- The effective 
theory decay rate is obtained by matching the current J onto the SCETn current given in 
Eq. plus a quark operator: 

(7(g)pi...Pn|JMJT(P)) 
- Er«/^^M(7(?)Pi---PnlxV-'^'^pTr{i?!C'f''^^)(P+;/i)i?a^1T(P)) 

i 

i 

= (0|xV-'^'^pl^(^))(7(9)l^10)^?i 



X 



(pi...p„|Tr{i?>f(P+;/i)Pf }|0) + {p,...Pn\L,,.W^ ^i(^;/^) w^i^^,M 



.(23) 



In obtaining the second line we make use of the results in Eqs. ()18|) and p9p. and use the 
properties of SCETn to factor soft and collinear degrees of freedom. In the last line we 
changed to a nonrelativistic normalization for the T state. 
Next we define the light-cone wave functions 



(pi . . .pn\VTi[Bl5^^r^Bt]\Q) = M^->g^"(x), 

fl (24) 
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where states are relativistically normalized, and the discrete label uj is converted to a contin- 
uous one, X = tu/n-p, as explained in Ref. [5|. The wave functions are dimensionless. See 
Appendix^for the relation of these SCET light-cone wave functions to those conventionally 
defined in QCD. Then the coUinear matrix elements in brackets in Eq. (|23|) can be written 
as the convolution: 



M 



2-n 



(25) 



The dependence on the scale fi cancels between the long-distance matching coefficients and 
the wave function. We will elaborate on this point in a moment. First we expand both 
; /i) and (f)^J^{x; ^) in Gegenbauer polynomials: 





[x-.jj) = 








n odd 


a{ 


^X, /i) 








n odd 




^X, /i) 








n odd 




[x;iJ,) = 


n odd 



(26) 

Then the convolution becomes an infinite sum of products of Gegenbauer coefficients: 

dx{af{x;iJ,)(l)g"{x;iJ,) + a? (x; /i)0j'"(x; /i)) 

E (/i/l^f + f3/l<^if^)bl:\f^)) , (27) 



1 

-1 



n odd 



where 



An) 
J 5/2 



n(n + l)(n + 2)(n + 3) 



(„) _ (n + l)(n + 2) 

J 3/2 ~ 



(28) 



9(rz + 3/2) ' '^/^ n + 3/2 

We now return to the question of the scale. Here we pick yU ~ A which minimizes 
logarithms in the wave function; however, large logarithms of M/A then appear in the 
Wilson coefficients. These large logarithms are summed using the renormalization group 
equations in SCET. This calculation was carried out in Ref. and we only quote the 
results here. We find: 



J o?x(af(x;/i)0j'"(x;/i) + a?(x;/i)0f"(x;/i)) 



il{M) J2 



n odd 



(n) 
7V 



An) (n) 

h/2 iW 

An) ^(n) 
h/2 



as{M) 



2aV"'//3o 
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in) 



a^(/i) 



2a(_"'//3o- 



as{M) 



asiM) 



(29) 
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where 



/3o 



7?^ 



(n) 
Iqq 



1 31 



'99 



The quantities \^^^ which appear in the exponents in Eq. (j29|) are negative for any n > 1. 

Furthermore A'^'' < 0, while A^"* = 0. This property allows us to consider the asymptotic 
hmit M > A, where as{M) 0. Then 
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2nt 



7^^ + it^ ± A(") 



(n) (n) 

Tgg - %q 



1 



+ 47J?7J? 



(n + l)(n + 2) 

+ 3n + 4 
l~(n + l)(n + 2) ' 
+ 3n + 4 

3n/ 



1 :^i 

2 > - 

9 ^ 7 



i=2 



n(n + l)(?2 + 2)(?2 + 3) 
2 2 



n(n + 1] 



+ 



> + 2)(n + 3) 



1 
6 



n+l 



j=2 



3 ' 



(30) 



lim / dx Taf (x; /i) + a^(x; /i)0f'"(x; 

M-+00 L y 



16 Cf 



a{{M) 



(31) 



which defines a nonperturbative parameter 5^". However, for values of M around the T 
mass this is not a very good approximation, and for values around the J/ip mass a much 
better approximation is to assume no running at all. 



A. Two body decay: T 7/2 

Having taken care of the technical details we can now use the above results to study the 
two body radiative decay T 7-Fi. The decay rate is 



r(T ^ ^F,: 



SCETt 



:(T|<,alx-p'xVl^p|T) 



167r'"'"P 

X J dx{al{x;n)(j)g'{x;n) + a'i{x;iJ,)(pg'{x]ix)) 



(32) 



where the full expression for the term in brackets is given in Eq. (|29p. After factoring, the 
soft matrix element involving the heavy quark fields was further simplified using the vacuum 
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insertion approximation for the quarkonium sector, which holds up to corrections of order 
|lO|. Note that the operator above overlaps only with the A = ±1 helicities of the T. 
Then using the rotation symmetries of NRQCD ^] we can relate the non-relativistic matrix 
element above to those conventionally used: 

(T|<alx-p'xt_p^l^p|T) = ^(T|0(l,3Si)|T) . (33) 

For the final state meson Fi to have nonzero overlap with the operators in Eq. ()24j) it must 
be flavor singlet, parity even and charge conjugation even. One candidate with the correct 
quantum numbers is the /2(1270). Furthermore this decay has been measured both in T 
and J/ip radiative decay, which is why we consider it. An interesting point is that only the 
helicity A = component of the /2 contributes at the order to which we are working. To 
see this begin by considering the decomposition of the following gluon matrix element into 
all possible light-cone form-factors: 

{h\ Tr[BlBi]\0) = A{e{\))gf + S^ef (A), (34) 

where e"'^ is the symmetric-traceless polarization tensor of the f2- We give the explicit form 
in Appendix ^ There are only two form factors above since the matrix element must be 
decomposed into tensors that have non-zero perpendicular components. The only structures 
available are g"^^ and . For A = ±1, e'^{X = ±1) = 0, so this helicity component 
does not appear at this order. The coefficient A{e{X)) is a scalar function which can be 
constructed from Tr(e_L) and fianpe'^^ . Because the helicity-zero polarization tensor has the 
property that e^^{\ = 0) oc , and the helicity-two polarization tensor has Tr(e^) = and 
UaUpe'^^ = 0, we can fix the normalization of the coefficient y4(e(A)) so that the first term on 
the right-hand side of Eq. (j34|) parameterizes the A = contribution while the second term 
parameterizes the A = ±2 contributions. The helicity-zero piece is picked out by the aig'^ 
term in Eq. ()18|) . while the helicity-two piece is picked out by the 02 term. Thus at leading 
order in perturbation theory, the dominant decay should be to the helicity-zero component 
of the /2. 

The NRQCD matrix element in Eq. ()32|) can be expressed in terms of the leptonic decay 
width of the T. At leading order, 

r(T e+e-) = ^^(T|0(l,='5i)|T) , (35) 
and the decay rate for T 7/2 can be expressed as: 

r(T 7/2) = ^6^^-(^)' (i^/.)2r(T e+e-) . (36) 

We can repeat the same analysis for the decay rate r{J/ip 7/2) and form a ratio of 
branching fractions, which in the asymptotic limit is: 



BiT ^ 7/2) 



bb) 



ACf + sY i?(T^e+e-) ^^^^ 



where Mqq = 2mQ. Using = 4.1 - 4.4 GeV, = 1.15 - 1.35 GeV, B{T e+e") = 
(2.38 ± 0.11) X 10-2, and 5(J/^ ^ e+e") = (5.93 ± 0.10) x lO'^ Q, we predict the ratio 
of branching fractions to be in the range: 

5(T-7/2) 



5(J/^^7/2)J 



= 0.14-0.19. (38) 
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As was mentioned earlier the asymptotic limit is not particularly good for the T, and quite 
bad for the J/ip. As a consequence we consider the scenario where the resummation of 
logarithms is neglected. In this case the ratio of branching fractions lies in the range 

BiT 7/2) , , 

^^7. =0.18-0.23. 39 

B[J/ip 7/2) 

We can improve this approximation by keeping more terms in the resummed formula in 
Eq. ()29|1 . The dominant term is the part of the n = 1 term proportional to h^g\ij)^ and in 
this approximation: 



B{T ^ 7/2) 
B{J/^ ^ 7/2) 



bb) 



0.13-0.18 



7?)-7i'^(«.(/i)M(M,^)) 



„(i) 



B{J/^ 



(40) 

where ~ 1 GeV. The range of values has not changed much from Eq. (jHH|l : however, 
theoretical errors are reduced: corrections to Eq. P0|l from the and higher-order terms 
in Eq. ()29j) are estimated to be roughly 40%, while corrections to the infinite mass limit 
from higher order terms are estimated to be roughly 80%. Corrections to Eq. (j39|) are hard 
to estimate; however, the range of values obtained give a rough upper limit on the ratio of 
branching ratios. In addition there are theory errors from neglecting higher-order terms in 
the perturbative expansion, as well as in the expansions in v and rj. Our prediction can be 
compared to the measured value of 0.06 ± 0.03, using the measurements B(T 7/2) = 
(8 ± 4) X 10"^ and B{J/iIj 7/2) = (1.38 ± 0.14) x lO'^ Given the theoretical errors 
we can only conclude that our prediction does not disagree with data. 

Our predictions for the ratios of T and J/ip branching fractions to 7/2 are consistent 



with the results of Refs. ji^, 2^, 25[ at leading order in the twist expansion used therein. In 
particular, we reproduce the suppression of the helicities |A| = 1, 2 in the final state relative 
to A = 0. In contrast with Ref. ^], we extract the NRQCD color-singlet matrix elements 
from the leptonic decay widths of T and J/ ip instead of the decay widths to light hadrons, 
for which corrections from color-octet contributions must be taken into account for a reliable 
calculation [s^. The leptonic decay width, however, receives large corrections at NNLO in 
perturbation theory [35|,|36|. In either case, one hopes that the uncertainties are mitigated 
in taking the ratios of branching fractions. 

We can also compare the decay rates of J/ip and ip' to 7/2 predicted by Eq. (jH^ at the 
matching scale = M, where af(x; M) = and af (x; M) is a constant. Dependence on the 
integral of the wave function </'{^(x; M) cancels out in the ratio of branching fractions: 



B[J/^ ^ 7/2) _ BjJ/j, ^ e+e- 
B{i;' ^ 7/2) 



B{i)' 



while the measured value is 6.57 ± 1.42. We used B{'il)' 
B{i)' 7/2) = (2.1 ± 0.4) X 10^^ 13^' 



= 7.85 ±0.35, (41) 
e+e-) = (7.55 ±0.31) x lO^^ and 



B. Three body decay: T — > •jinr 

Next we consider a two pion final state in the kinematic region where the pions are 
coUinear to each other with large energy and small total invariant mass ~ A. In this 
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case we we have a three body final state where the two pions are colhnear. It is convenient 
to define the variables: 



{P1+P2) 



n-pi 



M- 



T 



(42) 



In terms of these variables the differential decay rate is 



dr 



dml^ dz 5127r3M2 



X 



(ix(a?(a;;/i)0^''(x;/i) + a?(x; /i)0g''(a;; /i)) 



1 



(43) 



to leading order in m1^/M^. The p rope rties of the meson pair light-cone wave function (p'^'^ 
have been investigated in Refs. jSTi . [38| , which interestingly find that in the region where 
A <^ m,r7r ^ Mr they are given by an integral over two single-particle wave functions. The 
ratio of the T and J/ ip rates to 77r7r in the kinematic region of low m^^ is numerically the 
ratio in Eq. ()4Up times an extra factor of ml/ml ~ 0.07 — 0.1, that is, 0.01 — 0.02. This 
suppression is due to the much larger total phase space available in T — > 77r7r relative to that 
in T — »• 7/2. No T — > •yiin events have yet been observed in the region < (1.0 GeV)^ 



V. CONCLUSIONS 

We have systematically analyzed the exclusive radiative decays of quarkonium to ener- 
getic light mesons within the framework of soft-coUinear effective theory and non-relativistic 
QCD to leading order in the effective theory power counting, as well as to leading order 
in the strong coupling. We show that color-octet contributions are suppressed by a factor 
of f^as(/ic)/ ^/cis{M) ^ 0.05 in exclusive T decays, and can therefore be safely neglected. 
This is different from the situation in inclusive radiative decays in the endpoint region where 
octet contributions must be kept. 

We then turn to the color-singlet contribution. The tree-level matching onto this oper- 
ator is carried out in Refs. 0, 13 ! however, the authors do not consider the complete set of 
operators that could appear in this decay. We use the symmetries of SCET and NRQCD, 
including RPI, to show that the operator which is matched onto in Refs. [2, 3] is the only op- 
erator that can appear for the decays in question. We also consider the set of possible quark 
operators which can arise. Again only one of the possible quark operators can contribute 
to the decays we are interested in. This operator has zero matching coefficient, but it can 
be generated through running. We use the results of Ref. Q for the renormalization group 
mixing of the quark and gluon operators, thus resummin g la rge logarithms. Our results 
agree with the leading- twist analyses carried out in Refs. |23l 124 . Esj l . In Ref. jj25| higher 
twist corrections to the color-singlet contribution were considered. These corrections would 
be part of higher-order SCET corrections, which we have not studied here. Such an analysis 
would be complicated by the subleading color-octet contribution discussed in this work. 

Finally we study the phenomenology of quarkonium radiative decay to the /2, as well 
as to Tin where the pions are colhnear. We make predictions for the ratios of branching 
fractions B{T 7/2)/5(J/^ ^ 7/2) and B{J/^ lf2)/B{iP' 7/2), as well as for 
the differential decay rates of T and J /if) to 77r7r in the kinematic region of two colhnear 
pions. Our predictions for the decays to 7/2 are consistent with experimental data, but 
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with large theoretical uncertainties, while there is insufficient data for 77r7r with which to 
compare. Further theoretical work and more experimental data, especially for the light-cone 
wave functions of /2, will improve the precision of these predictions greatly. 
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APPENDIX A: NONPERTURBATIVE MATRIX ELEMENTS AND LIGHT- 
CONE WAVE FUNCTIONS 

The matrix elements in Eq. (j24|) defining the SCET wave functions can be related to 
conventional QCD wave functions for flavor- singlet mesons. The two-gluon wave functions 
for a meson with momentum q and net helicity A = 0, ±2 are defined as [ji^ ]: 

(0| TTG,,iz)Y{z, -z)G,,{-z)\q, A = 0),„ = f^q.qx £ e*^-^^0i(C, /io) , 
(0| Ti G,,{z)Y{z, -z)G,x{-z)\q, A = ±2)^„ 



where Y{z, ~z) is the path-ordered exponential of gluon fields: 



Y{z, —z) = Pexp 



ig J da ■ A(a) 



(Al) 



(A2) 



Going to the light-cone frame where = ^n^^ and z^^ = ^n^, we invert these formulas to 
find 



05(C;/^o) 



^M^A /-oo 



/oo 
e-C.-^+/2^0| Tr -;.+)G^(-z+)|g, A = 0) , 

-oo 

/oo 
-oo 

(A3) 



S ^ J -oo 
poo 



where z~^ = n-z, and q = n-q. Now we match the QCD fields on the right-hand side to 
fields in SCET: 



n^G^,{z+) 



Y(z^,-z^ 



(A4) 
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where 

Gr = -gP'- 'sA^n,,. - ^gAl,]. (A5) 

with 

= _p + + y2n-D„,. (A6) 

Therefore, for example, the matching between the QCD hght-cone wave-function 0^ and 
the SCET operator is 



'^i(C;/^o) - / rfz+e-<''-^"/2(0|TrG;:,(0)iy„e^^+^"/2iyiG"^(0)|g,A = 0) 

= _ E / rf^+e-^(^''--)^"/2(g--a;)(g- + ^)(0|Tr515p^,^5,^|g,A = 0). (AT) 

luTTJ^ 9 J — oo 



Integrating over and converting from the discrete index to a continuous where 
= LOcjq" we obtain the matching relation between the QCD and SCET light-cone wave 
functions 

05'^(C;/^) - -^(1 - C^)0f "'^HC;/^) . (A8) 

The SCET wave functions on the right-hand side are given by [cf. Eq. ()24|l ]: 

(0|PTr[5^5^,^^5;^]|M(g)) = r)^0f , 



(A9) 



Relations between the wave functions for the quark operator in QCD and SCET can be 
derived as in Ref. 130|. 



APPENDIX B: SPIN-2 POLARIZATION TENSORS 



The spin-2 polarization tensor for a particle of mass m and three-momentum k can be 
built from spin-1 polarization vectors using Clebsch-Gordan coefficients to arrive at 



e'^^(A = ±2) 



(A = ±i: 



/O o\ 

1 ±i 

±i -1 

yo 0/ 

/ Ik 



1 
2m 



|k| 
±i|k| 

/ k^ 



^^ = ^) = —.\ -o 



-in^ 

" 2 



±i\k\ \ 


ii^k 
±iEk J 

|k| Ek\ 






V|k|^k 














El 



J 



(Bl) 
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